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ABSTRACT 

We have investigated the effects of incorporating 
representations of heterogeneous chemical processes 
associated with stratospheric sulfuric acid aerosol 
into the LLNL two-dimensional, zonally averaged, 
model of the troposphere and stratosphere. Using 
distributions of aerosol surface area and volume 
density derived from SAGE II satellite observations, 
we were primarily interested in changes in 
partitioning within the Cl- and N- families in the 
lower stratosphere, compared to a model including 
only gas phase photochemical reactions. We have 
considered the heterogeneous hydrolysis reactions 

N 2 0 5 + II 2 0( aerosol) -> 2 IIN0 3 (T) 
and 

C10N0 2 + 11 2 0( aerosol) -4 IIOCl(T) + MN0 3 (t) 

alone and in combination with the proposed 
formation of nitrosyl sulfuric acid (NSA) in the 
aerosol and its reaction with IICl 

IIC1 + N0IIS0 4 ( aerosol) -> CINO(T) + II 2 S0 4 . 

Inclusion of these processes produces significant 
changes in partitioning in the NOy and ClOy families 
in the middle stratosphere. 

INTRODUCTION 

Numerical models of the troposphere and 
stratosphere including only gas-phase homogeneous 
photochemical processes can reproduce many of die 
observed features of trace constituent distributions. 
However, homogeneous models can not represent the 
chemical processes (Solomon, 1988) responsible for 
the observed lower stratosphere winter polar behavior 
of ozone and partitioning among the nitrogen and 
chlorine species. These processes require local 
temperatures low enough to condense water or a 
mixture of nitric acid and water into particles 
providing a surface for heterogeneous reactions 
which transform more stable chlorine reservoir 
species into more labile species. 

Homogeneous models also fail to produce the bias 
toward nitric acid observed in the partitioning of the 


nitrogen oxide family at mid and high latitudes in the 
winter hemisphere. Two reactions of importance that 
would take place on the surface of stratospheric 
sulfuric acid aerosol have been identified through 
laboratory investigation (Hofmann and Solomon, 
1989, and WMO/UNEP, 1992) 

N 2 0 5 + lI 2 0(acrosol) -> 2 HN0 3 (t) (1) 

and 

C10N0 2 + 1 1 2 0( aerosol) -> HOCl(T) + HN0 3 (T). (2) 

Both reactions convert NO x to IINO3, supplementing 
the homogeneous tcrmolccular reaction of OH with 
N0 2 . 

In consideration of the upcoming availability of a 
vastly increased data base of stratospheric trace 
species from the Upper Atmospheric Research 
Satellite (UARS), we have investigated the effects, in 
a numerical atmospheric model, that the sulfuric acid 
aerosol reactions above would have on species 
distributions, compared to inclusion of only gas phase 
processes. In this model study we have used an 
aerosol burden representing a climatological average 
of several years of SAGE II satellite data 
(WMO/UNEP, 1992), but, in order to distinguish the 
effects of the sulfuric acid aerosol reactions, have not 
included processes related to polar stratospheric 
clouds (PSC). Future direct comparison to UARS 
observations will be done with inclusion of the 
observed aerosol burden as perturbed by the Mt. 
Pinalubo eruption, as well as with PSC processes. 

We have also taken a preliminary look at the 
effects of an additional process proposed to occur in 
sulfuric acid aerosol of appropriate composition. 
Nitrosyl sulfuric acid (NOHSO4) has been tentatively 
identified as a common constituent of stratospheric 
sulfuric acid aerosol (Farlow et al., 1977) and is 
known from bulk phase laboratory work to be 
capable of existence as an ionic solute or solid within 
a range of sulfuric acid concentration and temperature 
characteristic of stratospheric aerosol. Burley and 
Johnston (1992) have shown that, based on current 
thermodynamic knowledge of the system, it is 
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possible that gaseous IIC1 in the stratosphere could 
react with NSA present in aerosol to produce nilrosyl 
chloride 

I1C1 + NOI IS0 4 (acrosol) -4 ClNO(t) + II 2 S0 4 . (3) 

CINO absorbs in die near ultraviolet with breaking of 
die Cl-N bond. Reaction (3), if it occurs, would dien 
act to shift die partitioning of C10 y toward Cl- 
conlaining radical species and temporary Cl 
reservoirs, such as C10N0 2 , and away from MCI. Its 
importance depends on its magnitude compared to the 
reaction of Oil widi MCI, which also releases Cl. 

MODEL DESCRIPTION 

The model used in diis study is a version of die 
Lawrence Livermore National Laboratory (LLNL) 
two dimensional model of the troposphere and 
stratosphere, with about 10° latitude rcsoludon and 3 
km vertical resolution in die stratosphere extending 
from pole to pole and die surface to 60 km. The gas 
phase homogeneous photochemistry in die model 
encompasses about 50 species and 120 diermal and 
photolylic reactions, for which kinetic and spectral 
parameters are taken from current recommendations 
(NASA, 1990). The circulation in the model is 
calculated interactively from the net healing rales 
resulting from die modeled species distributions and a 
climatological temperature field. In order to isolate 
the local photochemical effects of, for example, 
additional reactions, die circulation can be fixed to a 
previously calculated result. 

Reactions (1) and (2) arc treated as kinctically 
controlled, where the rate of the reaction is 
proportional to the rate of collision dirough the 
dimensionless reaction probability y. The first order 
rate constant for reaction is given by 



where v is die average molecular velocity and £ is 
die aerosol surface area density. In diis study, v is 
taken to be 5200 cm s' 1 for bodi N 2 0 5 and C10N0 2 
throughout die model domain. Based on laboratory 
results, die reaction probability for N 2 05 is assigned 
0.1. TorC10N0 2 the reaction probability is given by 

7(2) = 0.006cxp(-0.15(r-200)), 

where die temperature dependence arises from die 
dependence of die reaction probability on aerosol 
composition, which is in turn dependent on 
temperature. 

Surface area density for this study is based on an 
analysis of SAGE II data by Poole ct al. 
(WMO/UNEP, 1992) and is representative of an 
average aunosphcrc unperturbed by any major recent 
volcanic eruption. The prescribed values arc 
functions of altitude, latitude, and, coarsely, season, 
extending from 12 to 32 km and pole to pole. 


Minimum and maximum values are about 0.025 and 
1.75 * 10 ' 8 cm 2 cnr 3 , respectively. 

The rate constant for the proposed reaction of MCI 
with NS A dissolved or incorporated in the aerosol 
should depend on the NSA concentration in the 
aerosol, and may be controlled either by the 
collisional rate constant for MCI with the particle and 
an associated reaction probability, or by the Henry’s 
Law solubility for HC1 and a solution phase rate 
constant The presence of NSA in the aerosol depends 
on the concentration of H 2 S0 4 . According to Burley 
and Johnston (1992), NSA can be present at lower 
stratospheric temperatures when the H 2 S 04 wci ght 
per cent exceeds about 60%. The aerosol composition 
in the model was determined by interpolating from 
the tabulated results of Jaccker-Voirol (Jaecker- 
Voirol ct al., 1990) for binary H 2 O/H 2 SO 4 mixtures, 
using the model temperature and water partial 
pressure. Fig. 1 below shows the calculated values for 
July. The modeled concentration is generally below 
60% in the lowest portion of the aerosol distribution, 
but above 60% in the 15-30 km region of greatest 
importance to ozone controlling processes. 



Figure 1. Derived composition of July stratospheric sulfuric acid 
aerosol in H 2 SC >4 weight per cent. Solid contours are greater than 
60 per cent. Dashed contours are smaller than 60 per cent. 


Burley and Johnston (1992) propose several 
possible formation reactions for NSA, involving 
various NO y species, including 

NO + IINO 3 + H 2 S0 4 -4 NOIIS0 4 + N0 2 + II 2 0 (4) 
and 

NO 4 N0 2 + 2 II 2 S0 4 -4 2 N0IIS0 4 4 ll 2 o. (5) 

Farlow ct al (1977) tentatively observed solid NSA 
crystals in the collected aerosol, which would be 
formed when die NSA concentration in the aerosol 
exceeded its solubility in II2SO4. A proper 
calculation of the NSA abundance would require 
representing the processes of aerosol formation, 
growth, and loss, as well as the heterogeneous NSA- 
forming reactions. In this preliminary study we have 
assumed that the processes forming and destroying 
NSA arc rapid compared to its abundance, allowing 
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the use of a steady state expression for NSA 
abundance. The formation rate is taken as equal to the 
NO/aerosol collision rate times a reaction probability 
of 0.01, following Burley and Johnston (1992). The 
NS A loss rate is assumed to be the HCI/aerosol 
collision rale limes a reaction probability representing 
die relative abundance of NSA to II2SO4 at the 
particle surface, so that die atmospheric concentration 
of NS A is 


[NSA(jj)] 


f 

2*0.01*[NO]* 

V 


V^lhSOj} 

llici] 



V is die aerosol volume density and (U2SO4) is die 
concentration of sulfuric acid in the aerosol. V is 
simply sealed from die surface area density assuming 
dial die surface area in 1 cm 3 of aunosphcrc is on one 
spherical particle. Fig. 2 shows the resulting 
atmospheric concentration of NSA when these 
assumptions arc made. The values arc typically a few 
per cent or less of die H 2 S0 4 concentration. 
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Figure X Calculated rales of reaction in December (molecules 
cm' 3 s' 1 ) for (a) HCI + OH -4 Cl + Ii 2 0, ami (b) 1IC1 + NSA -> 
CINO + I^SC^. 

The model was integrated to steady state using 
1990 ambient boundary conditions for die source 
species, establishing the circulation for die case 
including reactions (l)-(3). This circulation was then 
prescribed for integrations to steady state for die case 
including reactions (1) and (2) and for the case of 
only gas phase reactions. 

RESULTS 

An expected major effect of both reactions (1) and 
(2) is a large increase in IINO3 (Fig. 4) at middle and 
upper latitudes particularly in winter and spring, 
where die rates of (l)+(2) become larger than die rate 
of Oil + N0 2 + M. This is shown in Fig. 5, where in 
December, the heterogeneous processes dominate 
nordiward of about 30°. 


Figure 2. Calculaled steady slate concent ration of NSA in mol 
cm* 3 , (a) December, (b) July. 


Given diese assumptions, the rate of I ICl loss and 
C1NO produedon in reaction (3) is expressed as 


djClNOj 

dl 


v S 

4 (V{h 2 SO.})^ 


[HCl|NSA(«)]. 


The calculated rate of MCI + NSA is compared to die 
rate of Oil + MCI in Fig. 3. 


As IINO3 increases, other NO y family members, 
particularly NO, N0 2 , and N2O5, show a 
corresponding decrease. The abundance of C10N0 2 
changes little because the y for reaction (2) is fairly 
small except in die soudicm high latitude winter. The 
effect on ozone of reactions (1) and (2) is to increase 
the importance of die C10 x destruction catalytic 
cycles while somewhat diminishing the dominant 
position of die NO x cycle between 20 and 30 km. 
The net overall effect on ozone is a slight increase of 
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0.1 1% relative to the total abundance in die gas phase 
ease. 



Figure 4. Derived IlN0 3 mole fraction at 62° N in December. 
Soliil line includes reactions (I) and (2), dashed line includes (1)- 
(3), dot-dash line includes only gas phase photochemistry. 

The addition of reaction (3) and die associated NSA 
formation and C1NO photolysis reactions do not 
produce significant further changes in NO y 
partitioning (Pig, 4). However, as expected, CIO and 




Figure 5. Calculated rates of HN0 3 production in December (mol 
cm ' 3 s' 1 ), (a) N0 2 + OH + M — ► HNO 3 + M, (b) Sum of N 2 0 5 + 
IWO(aerosol) -> 2 I!N0 3 and G0N0 2 + I l 2 0 (nerosol) -> HOG 
4-HN0 3 . 

C10N0 2 arc substantially increased between 20 and 
30 km globally, while I1C1 decreases. Pig. 6 shows 
dial in some regions of die stratosphere, widi the 


inclusion of reaction (3), C10N0 2 becomes the 
dominant inorganic Cl-containing species. Total 
ozone abundance is decreased by about 3% relative to 
the gas phase ease, widi the largest effects in regions 
where die photochemical ozone lifetime is long and 
the enhancement of CIO produces proportionately a 
larger reduction in local ozone abundance. 



Figure 6 . Calculated ratio of G0N0 2 to HG in December with 
NSA included. Where the value exceeds 1, C10N0 2 is the 
dominant inorganic G species. 


CONCLUSION 

This preliminary study of die effects, on trace 
species distributions, of heterogeneous processes 
taking place on stratospheric sulfuric acid aerosol 
provides some guidance for consideration of 
observations of NO y and C10 y family members. A 
strong enhancement of IINO3 in die partitioning in 
die NO y family at winter high ladludes is an indicator 
of aerosol moderated hydrolysis of N2O5, while NO, 
N0 2 , and N 2 0 5 arc substantially decreased. With die 
proposed NSA-rclatcd chemistry, CIONO2 is made 
die major inorganic Cl-containing species in die 24- 
31 km region. 

Puture improvements in our knowledge of the 
dicrmodynamic and kinetic properties of the sulfuric 
acid aerosol/H 2 0/NO y system will allow better 
model representation of NSA and related chemistry 
dian was attempted here. 

Consideration of heterogeneous chemistry on or 
within sulfuric acid aerosol is important for 
understanding of natural and anthropogenic 
influences on the stratosphere. Comparisons to 
observations of the UARS satellite and other data 
may provide evidence for die importance of die sc 
processes. 
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